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ABSTRACT 

Context. The physical structure of proto-planetary disks is not yet well constrained by current observations. Millimeter interferometry is an 



in 



essential tool to investigate young disks. 

Aims. We study from an observational perspective the vertical and radial temperature distribution in a few well-known disks. The surface 
density distribution of CO and HCO + and the scale-height are also investigated. 
^ ■ Methods. We report CO observations at sub-arcsecond resolution with the IRAM array of the disks surrounding MWC 480, LkCa 15 and 
' DM Tau, and simultaneous measurements of HCO + J=l— >0. To derive the disk properties, we fit to the data a standard disk model in which 
f"*) , all parameters are power laws of the distance to the star. Possible biases associated to the method are detailed and explained. We compare the 
properties of the observed disks with similar objects. 

Results. We find evidence for vertical temperature gradient in the disks of MWC 480 and DM Tau, as in AB Aur, but not in LkCa 15. The disks 
temperature increase with stellar effective temperature. Except for AB Aur, the bulk of the CO gas is at temperatures smaller than 17 K, below 
O the condensation temperature on grains. We find the scale height of the CO distribution to be larger (by 50%) than the expected hydrostatic 
scale height. The total amount of CO and the isotopologue ratio depends globally on the star. The more UV luminous objects appears to 
CZ3 . have more CO, but there is no simple dependency. The [ 13 CO]/[HCO + ] ratio is ~ 600, with substantial variations between sources, and with 



^ ■ 

Or 



radius. The temperature behaviour is consistent with expectations, but published chemical models have difficulty reproducing the observed CO 
quantities. Changes in the slope of the surface density distribution of CO, compared to the continuum emission, suggest of a more complex 
surface density distribution than usually assumed in models. Vertical mixing seems an important chemical agent, as well as photo-dissociation 
by the ambient UV radiation at the disk outer edge. 

Key words. Stars: circumstellar matter - planetary systems: protoplanetary disks - individual: LkCa 15, MWC 480, DM Tau - Radio-lines: 
stars 



1. Introduction 

CO rotation line observations of low-mass and intermediate- 
mass Pre-Main-Sequence (PMS) stars in Taurus-Auriga (~ 
140 pc, Kenyon et al., 1994) provide strong evidences 
that both TTauri and Herbig Ae stars are surrounded by 
large (R m « ~ 200 - 800 AU) Keplerian disk s (GM Aur: 
Koerner et ail d 19931). G G Tau: iDutrev etail i 1 994. MWC 480: 
Mannings et alj dl997l) ). Differences between TTauri and 
Herbig Ae d isks are observed in the very inner disk, very close 



to the star dDullemond et alj 1200 U [Monnier & Millan-G abet 



2002). However, detailed comparisons, namely based on re- 
solved observations instead of SED analysis, between Herbig 



Send offprint requests to: V.Pietu, e-mail: pietu@iram. fr 

1 Based on observations carried out with the IRAM Plateau de Bure 
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Ae disks and TTauri disks are too few to allow conclusions 
about the whole disk. In particular, there is no specific study on 
the influence of the spectral type on the disk properties. This 
can be achieved by comparing the properties of outer disks im- 
aged by millimeter arrays. 

The first rotational lines of CO (J=l-0 and J=2-l) mostly 
trace the cold gas located in the outer disk (R > 30 AU) since 
current millimeter arrays are sensitivity limited and do not al- 
low the detection of the warmer gas in the inner disk. In the 
last years, several detailed studies of CO interferometric maps 
have provided the first quantitative constrains on the physical 
propertie s of the outer disks. Among them, the method devel- 
oped by iDartois et alj 020031) is so far the only way to esti- 
mate the overall gas disk structure. IDartois et alj 020031) de- 
duced the vertical temperature gradient of the outer gas disk 
of DM Tau from a multi-transition, multi-isotope analysis of 
CO J=l-0 and J=2-l maps. The dust distribution and its tem- 
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perature can be traced by IR and optical observations and 
SEP modelling dP'Alessio et al.ll999tlDullemond & Dominik 
2004). However Pietu et al. ( 20061) recently demonstrated that 
SED analysis are limited by the dust opacity. Due to the high 
dust opacity in the optical and in the IR, only the disk surface 
is properly characterized. Images in the mm domain, where 
the dust is essentially optically thin, revealed in the case of 
MWC 480 a dust temperature which is significantly lower than 
those inferred from SED analysis and from CO images. With 
ALMA, the method described in this paper, which is today only 
relevant to the outer disk, will become applicable to the inner 
disk of young stars. 

In this paper, we focus on the gas properties deduced 
from CO data. The analysis of the millimeter continuum 
data, observed i n the meantime, is presented separately 
(iPietu et _ al.l 120061)). Here, w e generalize the method devel 



12 CO and in HCO + . H igh resolution 12 CO data was obtained 
by Simon et al. ( 2000b who used it to derive the stellar mass. 



oped by lDartois et a l. (2003) and we apply it to two Keplerian 
disks surrounding the TTauri star LkCa 15 (spec. type K5, 
Duvert et al.l2000l) and the Herbig Ae star MWC 480 (A4 spec- 
tral type, Simon et alj2000l) . By improving the method, we also 
discuss the estimate of disk scale heights. Then, we compare 
their properties with t hose surrounding Herbig Ae and TTauri 
disks such as AB AurdPietu etal.120 05^). HD 34282 dPietu et alj 
20031) . DM Tau dDartois et alj|2003l) . 

After presenting the sample of stars and the observations in 
Sect. 2, we describe the improved analysis method in Sect. 3. 
Results are presented in Sect. 4 and their implications discussed 
in Sect. 5. 



2. Sample of Stars and Observations 

2.1. Sample of Stars 

The stars were selected to sample the wider spectral type range 
with only a few objects. Table[T]gives the coordinates and spec- 
tral type of the observed stars, as well as of HD 34282 and 
GM Aur which were observed in 12 CO by IPietu et alJ J2003I) 



and iDutrey et all 0998) respectively. Our stars have spectral 
type ranging from Ml to AO. All of them have disks observed 
at several wavelengths and are located in a hole or at the edge 
of their parent CO cloud. 

DM Tau: DM Tau is a bona fide T Tauri star. 

Guilloteau & Dutreyl d 1994b first detected CO lines indi- 
cating Keplerian rota t ion us ing the IRAM 30 meter telescope. 
Guil loteau & Dutreyl d 1 998b have confirmed those findings by 
resolving the emission of the 12 CO J=l— >0 emission using 
IRAM Plateau de Bure interferometer (PdBI) and derived 
physical param eters of its protoplan etary disk and mass for 
the central star. Dartoi s et al. (2003) performed the first high 
resolution multi-transition, multi-isotope study of CO, and 
demonstrated the existence of a vertical temperature gradient 
in the disk. 



LkCa15: This is a C TT isola t ed fro m its parent cloud. 
Duvert et al] d2000b and lOi et all d2003l) observed the disk in 



MWC 480: The existence of a disk around this A 4 Herbig Ae 
star was first reported by Mannings et al. d 1997b . Using data 
from OVRO, they perform a modelling of t he gas disk consis- 
tent with Keplerian rotation. More recently. lsimon etaL (2000) 
confirmed that the rotation was Keplerian and measured the 
stellar mass from the CO pattern. 

ABAurigae: ABAur is considered as the proto-type of the 
Herbig Ae star (of spectral type AO) which is surrounded by 
a large CO and dust disk. The reflection nebula, obse rved i n 
the optical, ext e nds fa r away from the star ( Grady et al.ll 1999b . 
Semenov et alj d2005 ) used t he 30 -m to observe the chemistry 
of the envelope. Pietu et al. ( 20051) has analyzed the disk us- 
ing high angular resolution images of CO and isotopologues 
from the IRAM array and found that, surprisingly, the rotation 
is not Keplerian. A central hole was detected in the continuum 
images and in the more optically thin 13 CO and C ls O lines. 

2.2. PdBI data 

For DM Tau, LkCal5, and MWC 480, the 12 CO J=2-l data 
were observed in snapshot mode during winter 1997/1998 in 
D, C2 and B 1 configura tions. Detail s abou t the dat a quality and 



reduction are given in ISimon et al. I d2000h and in iDutrev et al 
1998). The 12 CO J=2-l data were smoothed to 0.2 km.s 



spectral resolution. The HCO + J=l— >0 data were obtained si- 
multaneously: a 20 MHz/5 12 channels correlator unit provided 



a spectral resolution of 0. 13 km.s 1 . The 13 CO J=l-0 and J=2 



1 observations of DM Tau are described in lDartois et a l. (2003). 

We obtained complementary 13 CO J=2-l and J=l-0 (to- 
gether with C ls O J=l-0 observations) of MWC 480 and 
LkCa 15. These observations were carried out in snapshot 
mode (partly sharing time with AB Aur) in winter periods be- 
tween 2002 and 2004. Configurations D, C2 and Bl were 
used at 5 antennas. MWC 480 has also shortly been observed 
in A configuration simultaneously with AB Aur. Finally, data 
was obtained in the largest A + configuration for LkCa 15 and 
MWC 48 in January 2006 (for a description of these observa- 



tions, see Pietu et al. 2006b . Both at 1.3 mm and 2.7 mm, the 



tuning was double-side-band (DSB). The backend was a cor- 
relator with one band of 20 MHz, 512 channels (spectral res- 
olution 0.11 km.s ) centered on the 13 CO J=l-0, one on the 
C ls O J=l-0 line, and a third one centered on the I3 CO J=2->1 
line (spectral resolution 0.05 km.s -1 ), and 2 bands of 160 MHz 
for each of the 1.3 mm and 2.7 mm continuum. The phase and 
secondary flux calibrators were 0415+379 and 0528+134. The 
rms phase noise was 8° to 25° and 15° to 50° at 3.4 mm and 
1 .3 mm, respectively (up to 70° on the 700 m baselines), which 
introduced position errors of < 0.05". The estimated seeing is 
about 0.2". 

Flux density calibration is a critical step in such projects. 
The kinetic temperature which is deduced from the modelling 
is directly proportional to the measured flux density. This is 
even more problematic since we compare data obtained on sev- 
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Table 1. Star properties 



Source 


Right Ascension 


Declination 


Spect.Type 


Effective Temp.(K) 


Stellar Lum.( L Q ) 


CO paper 


LkCa 15 


04:39:17.790 


22:21:03.34 


K5 


4350 


0.74 


1,2 


MWC480 


04:58:46.264 


29:50:36.86 


A4 


8460 


11.5 


1.2 


DMTau 


04:33:48.733 


18:10:09.89 


Ml 


3720 


0.25 


3 


AB Aur 


04:55:45.843 


30:33:04.21 


A0/A1 


10000 


52.5 


5 


HD 34282 


05:16:00.491 


-09:48:35.45 


A1/A0 


9440 


29 


4 


GM Aur 


04:55:10.98 


30:21:59.5 


K7 


4060 


0.74 


6 



Col. 2,3: J2000 coordinates deduced from the fit of the 1.3 mm continuum map of the PdBI. Errors bars on the as t romet ry ar e of order < 0.05" 
Col. 4,5, and 6, th e spectral type, effective te mperature and the stellar luminosity are those given in ISimon et alj feOOd) a ndlPietu et all J2003 
for HD 34282 and van den Ancker et al. I jl998h for AB Aur. Col . 6, list of C O interferometric pap ers are: 1 = this paper, 2 = ISimon et al J ( 12000 
3 = lDartois et at] ilOOA 4 = |Pietu et alj d2003n . 5 = |Pietu et al.U2005h . 6 = |Putrev et ail Jl998 



eral years. Hence, we took great care of the relative flux cali- 
bration by using not only MWC 349 as a primary flux refer- 
ence, assuming a flux S(v) = 0.95(v/87GHz) a60 (see the IRAM 
flux report 14) but also by having always our own internal ref- 
erence in the observations. For this purpose we applied two 
methods: i) we used MWC 480 as a reference because its con- 
tinuum emission is reasonably compact and bright, ii) we also 
checked the coherency of the fl ux calibration on the spectral in- 
dex of DM Tau, as described in lDartois et alj d2003l) . We cross- 
checked all methods and obtained by these ways a reliable rel- 
ative flux calibration from one frequency to another. 



2.3. Data Reduction and Results 



We used the GILDAS0 software package (Pety 2005) to re- 
duce the data and as a framework to implement our minimiza- 
tion technique. Figure Q] presents some channel maps of 12 CO, 
13 CO J=2-^l and 13 CO J=1^0 for LkCal5 and MWC 480. 



Simila r figures for DM Tau can be found in iDartois et al 
(2003), and for ABAur in IPietu et all d2005l) . A velocity 



smoothing to 0.6 km.s" 1 (for MWC 480) and 0.4 km.s -1 (for 
LkCa 15) was applied to produce the images displayed in Fig[TJ 
However, in the analysis done by minimization, we used spec- 
tral resolutions of 0.20 km.s" 1 ( 12 CO and HCO + ) and 0.21 
km.s -1 ( 13 CO), and natural weighting was retained. 

The HCO + images are presented in Fig 12 The typical an- 
gular resolution of these images is 3-4". For LkCa 15, the line 
flux is compatible w ith the result quoted bv lDuvert et al. I d2000h 
and lOi et"aT| d2003l). and for DM Tau with the single-dish mea- 
surement of Dut rey et al] dl9971) . 



3. Model description and Method of Analysis 

In this section, we summarize the model prope rties and de- 
scribe the improvements since the first version dDutrev et al 
1994. 



3.1. Solving the Radiative Transfer Equation 

To solve the radiative transfer equation coupled to the statisti 
cal equilibrium, we have written a new c ode, offering non LTE 
capabi lities. More details are given in IPavlyuchenkov et al 



1 See http://www.iram.fr/IRAMFR/GILDAS for more informa- 
tion about the GILD AS softwares. 



(2007). The LTE part, which is sufficient for low energy level 
of rot ational CO lines (J= l-0 and J=2-l), is strictly similar to 
that of lDutrev et al. (1994), and uses ray-tracing to compute the 
images. The HCO + data was also analyzed in LTE mode (see 
Sect |3.2.3l for the interpretation of this assumption). We use a 
scheme of nested grids with regular sampling to provide suf- 
ficient resolution in the inner parts of the disk while avoiding 
excessive computing time. 

3.2. Description of the physical parameters 

For each spectral line, the relevant physical quantities which 
control the line emission are assumed to vary as power law as 
function of radius (see Sect |3.2.1| for the interpretation of this 
term) and, except for the density, do not depend on height above 
the disk plane. The exponent is taken as positive if the quantity 
decreases with radius: 

a(r) = a (r/R a )- e ° 

For each molecular line, the disk is thus described by the fol- 
lowing parameters : 

- Xo, Yo, the star position, and Vdisk, the systemic velocity. 

- PA, the position angle of the disk axis, and i the inclination. 

- Vo, the rotation velocity at a reference radius R v , and v the 
exponent of the velocity law. With our convention, v = 0.5 
corresponds to Keplerian rotation. Furthermore, the disk is 
oriented so that the Vq is always positive. Accordingly, PA 
varies between and 360°, while i is constrained between 
-90° and 90° (see fig.[3>. 

- T m and q,„, the temperature value at a reference radius Rt 
and its exponent (see Sect l3.2.3l >. 

- dV, the local line width, and its exponent e,.. The interpre- 
tation of dV is discussed in Sect |3.2.6l 

- S m , the molecular surface density at a radius R% and its ex- 
ponent p m 

- R out , the outer radius of the emission, and R[ n , the inner 
radius. 
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Fig. 1. Channel maps of the 12 CO and 13 CO line emission towards MWC 480 and LkCa 15. The synthesized beam is indicated in 
each panel. The cross indicates the direction of the major and minor axis of each disk, as respectively derived from the analysis 
of the line data. The velocity of the channel is indicated in the upper left corner of each panel (LSR velocity in km.s 1 ). A 
200 m taper has been applied on the I3 CO J=2— >1 data of LkCa 15 which where obtained at higher angular resolution. Left: 
MWC 480 Left: 12 CO J=2^1 line, spatial resolution 1.38 x 0.99" at PA 16°, contour spacing 200 mJy/beam, or 3.4 K, or 3.6 
cr. Middle: 13 CO J=2->1 line, spatial resolution 0.77 x 0.57 at PA 26°, contour spacing 50 mJy/beam, or 2.9 K, or 1.9 cr. Right: 
13 CO J=l— >0 line, spatial resolution 1.72 x 1.24" at PA 37°, contour spacing 20 mJy/beam, or 0.9 K, or 1.8 cr. Negative contours 
are dashed, and the zero contour is omitted. Right: LkCa 15 Left: 12 CO J=2— »1 line, spatial resolution 2.01 x 1.20" at PA 12°, 
contour spacing 150 mJy/beam, or 1.4 K, or 2.5 cr. Middle: 13 CO J=2— >1 line, spatial resolution 1.46 x 1.24 at PA 52°, contour 
spacing 100 mJy/beam, or 1.4 K, or 1.6 cr. Right: 13 CO J=l— >0 line, spatial resolution 2.55 x 1.51" at PA 47°, contour spacing 
30 mJy/beam, or 0.8 K, or 1.8 cr. Negative contours are dashed, and the zero contour is omitted. 



- h m , the scale height of the molecular distribution at a ra- 
dius Ri u and its exponent e/,: it is assumed that the density 
distribution is Gaussian, with 

n{r, z ) = -^p exp [- (z/h(r)) 2 ] (1) 
h(r) \7i L J 

(note that with this definition, e/, < in realistic disks) 

thus giving a grand total of 17 parameters to describe the emis- 
sion. 



It is important to realize that all these parameters can ac- 
tually be constrained for the lines we have observed, under the 
above assumption of power laws. This comes from two specific 
properties of proto-planetary disks: i) the rapid decrease of the 
surface density with radius, and ii) the known kinematic pat- 
tern. The only exception is the inner radius, R[ n , for which an 
upper limit of typically 20 - 30 AU can only be obtained. 

The use of power laws is appropriate for the velocity, pro- 
vided the disk self-gravity remains small. Power laws have 
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Fig. 2. Channel maps of the HCO + J=l— »0 line emission towards DMTau, LkCa 15 and MWC480. The cross indicates the 
direction of the major and minor axis of each disk. The velocity of the channel is indicated in the upper left corner of each 
panel (LSR velocity in km.s 1 ). For displaying purpose, the spectral resolution is 0.6 km.s 1 for MWC 480 and LkCa 15, and 0.4 
km.s 1 for DMTau. All contour steps are 20 mJy/beam (~ 2<x), negative contours are dashed and the zero contour is omitted. 
Top: DM Tau, spatial resolution is 4.3 x 3.3" at PA 47°, contour step 0.22 K. Middle: LkCa 15, spatial resolution is 4.6 x 3.4" at 
PA 19°, contour step 0.25 K. Bottom: MWC 480, spatial resolution is 3.5 x 2.5" at PA 42°, contour step 0.36 K. 



been shown to be a good ap proximation for the kinetic tem- 
perature distribution (see, e.g. lChiang & Goldreichll 19971) . and 
thus to the scale height prescription. For the surface density, 
power laws are often justified for the mass distribution based 
on the a prescription of the viscosity with a constant accre- 
tion rate: self-similar solutions to th e disk evolution impl y 



power la ws with an ex ponential edge (Ha rtmann et alJ l!998). 



Although Mal bet et alJ (120011) pointed out some limitations of 
this assumption for the inner parts (r < 30 AU) of the disk, the 
approximation remains reasonable beyond 50 AU. However, 
although this may be valid for H2, chemical effects may lead 
to significant differences for any molecule. Given the limited 
spatial dynamic range and signal to noise provided by current 
(sub-)millimeter arrays, using more sophisticated prescriptions 
is not yet justified, and power laws offer a first order approxi- 
mation of the overall molecular distribution. 



and a cylindrical representation, where the functions depends 
on the projected distance on the disk (r = p) and height z 



F = f(p,z) 



(3) 



with d 2 = p 2 + z 2 ■ It is obvious that the exponent of the power 
laws will depend slightly on which description is selected. We 
have evaluated the differences between the two different repre- 
sentations (see Tab|2]): as expected, the results are very similar 
for 13 CO, but the cylindrical representation is a ~ 8 sigma bet- 
ter representation for 12 CO J=2— >1 (the most affected transi- 
tion because of its higher opacity). The most affected param- 
eters are S and h. However, with the exception of the scale 
height, the effects remain small compared to the uncertainties. 
Accordingly, we have arbitrarily decided to represent all pa- 
rameters in terms of cylindrical laws, i.e. the radius is r — p in 
the model description. 



3.2.1. Spherical vs Cylindrical representation 

To describe the disk and its parameters, we can use either cylin- 
drical or spherical coordinates. 

The description given in the previous sub-section is unam- 
biguous for thin disks (where h(r) <sc r), However, for molec- 
ular disks with outer radius as large as 800 AU, the disk thick- 
ness is no longer small, and the classical theory of hydrostatic 
scale height, which uses cylindrical coordinates, is only valid 
in the small h/r limit. It may then become significant to distin- 
guish between a spherical representation, where the functions 
depends on the distance to the star (r = d) and height above 
disk plane, z 

F = f(d,z) (2) 



3.2.2. Reference Radius 

Unless p m =i 0, the choice of the reference radius R% will affect 
the relative error bar on E m , (ffi,„/E m . Depending on the an- 
gular resolution and on the overall extent of the emission, for 
each molecular line, there is a different optimal radius Rz which 
minimizes this relative error. The same is true for all other ref- 
erence radii for the power laws. This makes direct comparison 
at an arbitrary radius not straightforward. To circumvent this 
problem, we have determined S m , &E m as a function of the ref- 
erence radius for each transition, and selected the R% giving the 
best S/N ratio. Results are given Fig|4]|6] note that the curves 
^m(Rz) should not be interpreted as independent estimates of 
the local surface density at different radii: they explicitly rely 
on the assumption of a single power law throughout the whole 
disk. 
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Table 2. Cylindrical vs Spherical description 





Fig. 3. Geometrical convention: four disk configurations yield- 
ing the same aspect ratio when projected are presented. The 
shaded area correspond to the part of the disk that is closer 
to us than the plane of the sky containing the disk axis. The 
curved arrow represents the sense of rotation (this defines the 
disk axis). The projected velocities allow to disentangle be- 
tween one of the two columns (for the two cases presented on 
the left, the red-shifted part would be the lower part, and the 
contrary for the two cases on the right). It is possible to deter- 
mine the sign of the inclination either with the asymmetry at the 
systemic velocity if sufficient signal to noise is available, or by 
other means, e.g. scattered images in the optical and near-IR. 
PA is classically designed to be positive East from North. 



The same arguments can be applied to all other power laws, 
in particular to the scale height and the temperature. However, 
for the latter, the exponent is usually small (q - -0.2 - 0.5), 
so the choice of the reference radius is less critical. Figure |4] 
shows the variation of the parameters (here S m , T m and li) and 
of the signal to noise ratio as a function of the reference ra- 
dius for DM Tau. Th is is a reanalysis of the data published by 
Dartois et al.l (120031) . This method gives consistent results with 
previous work, but allows a more precise determination of the 
parameters (especially the surface density as stated above), by 
determining in which region they are constrained. Figj4] indi- 
cates that the temperature is determined around 100 - 200 AU, 
while the surface densities are constrained in the 300 - 500 AU 
region. Somewhat smaller values apply to the other sources, 
LkCa 15 and MWC 480 (see Fig US). 



Line 


1J CO(1-0) 


13 CO(2-l) 


12 CO(2-l) 


Cylindrical 


Rout (AU) 


791 ±33 


725 ± 24 


868 ± 15 


r„(K) 


9.1 ±0.3 


14.7 ±0.3 


24.2 ± 0.2 


q 


-0.17 ±0.09 


0.18 ±0.05 


0.46 ±0.01 


Vo(km.s-') 


2.12 ±0.07 


2.02 ± 0.06 


2.09 ± 0.06 


V 


0.5 


0.5 


0.5 


dV Coils -1 ) 


0.133 ±0.015 


0.136 ±0.015 


0.147 ±0.010 


e v 


-0.28 ± 0.08 


-0.04 ± 0.07 


0.08 ± 0.04 


Z,„ (cm -2 ) 


26.9 ± 1.7 10 16 


41.4 ±6.1 10 16 


37.4 ±9.1 10 18 


Pm 


3.3 ±0.3 


3.7 ± 0.4 


4.5 ± 0.5 


h m (AU) 


27.1 ±4.4 


27.8 ± 2.9 


33.4 ±2.3 


eh 


-1.07 ±0.20 


-1.09 ±0.10 


-0.99 ± 0.04 


x 1 


513329.3 


514557.2 


136376.5 


Spherical 


Rom 


804 ± 34 


737 ± 20 


898 ± 12 


To 


9.3 ±0.3 


15.2 ±0.4 


26.7 ± 0.2 


1 


-0.14 ±0.09 


0.20 ± 0.05 


0.50 ± 0.01 


V (km.s- 1 ) 


2.16 ±0.07 


2.07 ± 0.06 


2.08 ± 0.05 


V 


0.5 


0.5 


0.5 


dV (km.s- 1 ) 


0.135 ±0.016 


0.137 ±0.015 


0.183 ±0.012 


e v 


-0.26 ± 0.08 


-0.05 ± 0.06 


0.17 ±0.05 


L m (CITT 2 ) 


25.7 ± 1.8 10 16 


35.1 ±5.3 10 16 


41.1 ± 7.4 10 17 


Pm 


3.2 ±0.3 


3.6 ±0.4 


3.4 ±0.4 


h m (AU) 


20.7 ± 2.2 


23.4 ± 2.4 


19.7 ± 1.4 




-1.32 ± 0.18 


-1.13 ±0.12 


-1.25 ±0.06 


x 1 


513329.5 


514558.9 


136451.5 



Results for DM Tau with "free" scale-height and "free" line-width (see 
Sect. 13.2.51 and [3"2"6l for more details) in the cylindrical and spherical 
representations. Caution: this table is here to illustrate the dependence 
on the geometry of as many parameters as possible. The best parame- 
ters to be applied to DM Tau should be taken from TablefJ] 



3.2.3. Temperature law 

The interpretation of T m (and as a consequence of £,„) depends 
on the specific model of radiative transfer being used. In this 
paper, we are dealing with CO isotopologues, for which LTE 
is good approximation: T m is then the kinetic temperature. For 
transitions which may not be thermalized, two different radia- 
tive transfer models may be used. We can apply an LTE ap- 
proximation: since we are fitting brightness distributions, T m is 
in this case the excitation temperature T ex , and the surface den- 
sity E m is computed assuming the same temperature controls all 
level population, i.e. the partition function. We can also solve 
for molecular line excitation using a non-LTE statistical equi- 
librium code: T m is then the kinetic temperature, and £,„ the 
total molecule surface density, within the limitations of the ra- 
diative transfer code accuracy. In this paper, only the LTE mode 
was used. 



3.2.4. Surface Density 

The surface density is expected to fall off rapidly with distance 
from the star — 1-2), and in particular much more 

rapidly than the temperature (q st 0.5). Remembering that, be- 
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cause of the partition function, the line opacity of a J — 1-0 
transition scales as Z, m /T 2 at high enough temperatures (for 
constant line width), this indicates that, in general, the central 
part of the disk is much more optically thick than the outer 
regions. For the detectable lines, the temperature can be de- 
rived from the emission from this optically thick core, while 
the surface density E m is derived from the optically thin region. 
This remains valid provided the same temperature law applies 
to the two regions. For the weakest lines, the optically thick 
core may be too small, and the temperature exponent remains 
essentially unconstrained: this should be reflected to by the er- 
ror bars on the temperature. In such cases, however, the line 
emission scales as E m /T (for J=l-0 line). 

If the molecule is very abundant, such as CO, the local 
line core may be optically thick almost throughout the disk. 
However, the local line wings remain optically thin: if the ve- 
locity resolution is high enough, this information can be used to 
constrain the opacity, and hence the surface density, although 
less accurately. To illustrate this point, let us consider the case 
of a face on disk. The apparent width of the line profile is then 
controlled only by the internal velocity dispersion (turbulent + 
thermal) and by the line opacity. In the outer parts, the line is 
optically thin and of Gaussian shape, while in the inner parts, 
the line is broadened by the large opacity, and the profile tends 
towards a square shape, with width proportional to y/ln(r). The 
evolution of the line profiles as function of radius indicates at 
which radius t = 1, and thus constrains the molecular sur- 
face density. Note that if the spectral resolution is insufficient 
to properly sample the line width, this is no longer possible. 
Current correlators are not instrumentally limited but smooth- 
ing can be necessary when the sensitivity is the limiting factor. 
If the disk is not seen face on, the systematic velocity gradi- 
ent get superimposed to this local line broadening, but does not 
change the fundamental relationship between the apparent (lo- 
cal) line width and the opacity. In any case, since the effect goes 
as s/ln(T), measuring the surface density by such a method is 
expected to be rather unprecise, and will result in larger errors. 
Furthermore, it should be noted that this opacity broadening 
effect results in a slight coupling of the line width parameters, 
dV and e v , to the surface density. 



3.2.5. Scale Height 

The existence of a known kinematic pattern also allows to con- 
strain the scale height parameters /zq and <?/,. Unless the disk is 
seen face on, for an optically thin line, the velocity range in- 
tersected along a line of sight will depend on the disk inclina- 
tion and flaring, resulting in a coupling between scale height 
and line width. For an optically thick line, the difference in 
inclinations between the front and back regions of the disk 
can also be measured. The spatial distribution of the emission 
as a function of projected velocity (or function of the veloc- 
ity channels) along the line of sight will depend on the scale 
heig ht of the disk (this is illustrated for the opacity by Fig. 3 
from iDartois et al. (2003)). At first glance, the expected scale 
heights, of order 50 AU at 300 AU from the star, would ap- 
pear too small compared to the resolution of 0.7 - 1.5" (100 to 



200 AU), and the exponent e/, even more difficult to constrain. 
However, millimeter interferometry, being an heterodyne tech- 
nique, allows phase referencing between velocity channels, and 
the precision in relative positions is equal to the spatial resolu- 
tion divided by the Signal to Noise (as our bandpass accuracy, 
about a degree, is not the limiting factor), i.e. ~ 10 AU. This 
super-resolution approximately matches the expected displace- 
ments due to Keplerian rotation, 



5r = 2r— 10 AU at r = 100 AU 



(4) 



for 5v — 0. 15 km.s -1 . Note that in Eq|U the value to be used for 
6v is approximately the largest of the local line width and the 
spectral resolution. As a result, the scale height has a measur- 
able effect on the images, even though the flaring parameter is 
difficult to constrain. 

Actually, because of the effect of the scale height on the 
disk images, if a simple hydrostatic equilibrium is assumed in- 
stead of handling the scale height separately, it will result in a 
bias in the temperature law since one is then introducing a cou- 
pling between the vertical distribution of the molecules and the 
temperature which is vertically isothermal in our model. 

Finally, note that because the derivation of the scale height 
depends on the differing inclinations of the optically thick 
parts, there is also a coupling between the surface density and 
the scale height derivation. This coupling is specially strong for 
disks close to edge-on. 

3.2.6. Line Width 

Line broadening results from a combination of thermal and tur- 
bulent velocities. 



AV 



2kT 



+ v 



turb 



(5) 



However, an accurate representation of the thermal broaden- 
ing requires a correct fitting of the gas temperature. This can 
become problematic when dealing with lines which are sub- 
thermally excited. Accordingly, our code offers two different 
parameterizations of the line width. The first one is 



AV 



2kT 



+ dV 2 i.e. dV = v t urb 



(6) 



which corresponds, for e Y = ( Vturb = cte), to the parametriza- 
tion used in our previous papers. The second one is 



AV - dV - dV m (r/R dv y e * 



(7) 



i.e., we do not separate the thermal component from the total 
line width. This offers the advantage of providing values which 
are not biased when the kinetic temperature is difficult to con- 
strain (as can happen for lines which are not at LTE). We use 
this new description in this paper. 

Note that all line widths we quote are half-width at 
1 /e. This allows direct comparison with the thermal velocity, 
y/2kT/m, but one should multiply by 2 Vln(2) = 1 .66 to convert 
to FWHM. As described before, the line width is also coupled 
to the scale height. 
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Fig. 4. From left to right 12 CO J=2->1, 13 CO J=2->1, 13 CO J=l-»0 and HCO + J=1^0 in DMTau. Top: Temperature T m 
(black curve, left axis) and signal to noise on the temperature T m /6T m (grey curve, right axis) as a function of reference radius 
Rj .Middle Surface density ~L m (left axis) and signal to noise on the surface density E m /6E m (right axis) as a function of reference 
radius R%. Bottom Scale height h (left axis) and signal to noise on the scale height h/Sh (right axis) as a function of reference 
radius /?/,. 



3.3. Minimization Technique 

Using the disk parametrization described above, we have im- 
proved several aspe cts of the method originally developed by 



Dartois et al. (2003). 



3.3.1. Minimization Method and Error Bars 

The comparison between model and data and the x 1 mini- 
mization is alw ays performed inside the UV plane, using nat- 
ural weighting (Gui lloteau & Dutrey Il998h . We have imple- 
mented a modified Levenberg-Marquardt minimization scheme 
to search for the minimum, and we use the Hessian to derive 
the error bars. This scheme is much faster than the grid search 
technique used previously, and less prone to local minima. It 
allows to fit simultaneously more parameters, and thus provide 
a better estimates of the errors because the coupling between 
parameters is taken into account. 

However, the current fit method does not handle asymmet- 
ric error bars which happen in skewed distributions: this limi- 
tation should not be ignored when considering the errors on S, 
T or R out . The quoted error bars include thermal noise only, but 



do not include calibration errors (in phase and in amplitude). 
Amplitude calibration errors will directly affect the absolute 
values of the temperature or the surface density (as a scaling 
factor, but not the exponent) while phase errors can introduce 
a seeing effect. The latter effect is very small for line analysis 
because molecular disks are large enough but it starts to be sig- 
nificant for dust disks. Our careful calibration procedure brings 
the amplitude effect to below 10 %. All other parameters are 
essentially unaffected by calibration errors. 

3.3.2. Continuum Handling 

The speed improvement allowed us to treat more properly the 
continuum emission. Even though the continuum brightness 
is small compared to the line brightness, continuum emission 
cannot be ignored when trying to retrieve parameters from the 
line emission because it appears systematically in all spectral 
channels. 

A simultaneous fit of the continuum emission with its own 
physical parameters will properly take into account these bi- 
ases. We performed this by adding to the spectral UV table a 
channel dedicated to the continuum emission. The global x 2 is 
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Fig.5. From left to right 12 CO, 13 CO and HCO + J=l->0 in LkCa 15. Top: Temperature T m (black curve, left axis) and signal 
to noise on the temperature T m /ST m (grey curve, right axis) as a function of reference radius Rj- Middle Surface density £ m (left 
axis) and signal to noise on the surface density H m /SE m (right axis) as a function of reference radius R^. Bottom Scale height h 
(left axis) and signal to noise on the scale height h/dh (right axis) as a function of reference radius R/,. For HCO + the temperature 
and surface density were adjusted separately: the curves should be used as an indicator of the region over which these values are 
actually constrained, but not as a quantitative measure in terms of S/N, as the coupling between temperature and surface density 
is ignored. 



given by the sum x 2 — X% + x\ wnere X% anc ^ x\ are tne X 2 f° r 
the continuum and for the line minimization, respectively. The 
continuum emission is fitted using the following parameters: 
Ram, Remit, the inner and outer radii of the dust distribution, 2^ 
and pd, the surface density power law, and 7^ and qj, the dust 
temperature power law. Note that the representation of the con- 
tinuum needs only be adequate within the noise provided by the 
spectral (total) line width, about 10 to 20 MHz, rather than for 
the whole receiver bandwidth (600 MHz). A simplified model 
is thus often acceptable, as shown below. 

We also compared this robust method with a simplified (and 
faster) one in which we subtract, inside the UV plane, the con- 
tinuum emission from the line emission before fitting. This sub- 
traction is not perfect: for the optically thick regions (which 
always occupy some fraction of the line emission) continuum 
should not be subtracted. However, unless the disk is seen face 
on, because of the Keplerian shear, at each velocity, the line 
emitting/absorbing region only occupies a small fraction of the 
disk area. This fraction is at most of order of the ratio of the 
local line width to the projected rotation velocity at the disk 



edge, in practice less than 10 - 15 % for the disks we con- 
sidered. Accordingly, the error made by subtracting the whole 
continuum emission in every channel is quite small. 

We verified that this simplified method gives, as expected, 
the same results as the robust one. We used it, as it is sig- 
nificantly faster. This better handling of the continuum al- 
lowed us to impro ve on the mass determination of MWC 480 
dSimon et alB oOO), where the continuum is relatively strong. 

3.3.3. Multi-Line fitting 

The code has also been adapted to allow minimization of more 
than one data set at a time. This has been used to simultane- 
ously fit the 13 CO J=l— >0 and J=2 — > 1 transitions. For multi- 
line fitting, the temperature is derived from the line ratio. The 
temperature is determined in the disk region where both lines 
are optically thin. In this case, the fit is much more accurate 
than from a single line observation if there is no significant 
vertical temperature gradient. When a temperature gradient is 
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Fig. 6. From left to right 12 CO J=2->1, I3 CO J=2->1, I3 CO J=l->0 and HCO + J=l-»0 in MWC480. Top: Temperature T m 
(black curve, left axis) and signal to noise on the temperature T m /6T m (grey curve, right axis) as a function of reference radius 
R j. Middle Surface density E m (left axis) and signal to noise on the surface density E m /6E m (right axis) as a function of reference 
radius R%. Bottom Scale height h (left axis) and signal to noise on the scale height h/dh (right axis) as a function of reference 
radius /?/,. As for LkCa 15, temperature and surface density are adjusted separately for HCO + . 



present, the derived temperature will reflect an "average" tem- 
perature weighted by the thermal noise. 

4. Results 

4.1. Interpretation of the Power Law model 

The parametric model would perfectly describe the line emis- 
sion from a molecule in LTE in a vertically isothermal disk 
which is in hydrostatic equilibrium, with power law for the ki- 
netic temperature and surface density, and constant molecular 
abundance. In such a case, e v = q/2 if the local line width is 
thermal, and h - q/2 - 1 - v and ho = ( ^2kTolm)IVo (us- 
ing the same reference radius Ri, = Rt = R v , m being the 
mean molecular weight). If such a disk was chemically ho- 
mogeneous, all molecular lines would yield the same results 
for the 17 parameters (after correction of £,„ for the molecule 
abundance). Differences between the parameters derived from 
several transitions will reflect departures from such an ideal sit- 
uation. 

For example, the geometric parameters PA and i should all 
be identical, as should be the kinematic parameters Vdisk, Vo, v 
(and we should have v = 0.50 for Keplerian rotation). Xq, Yq 



should reflect the absolute astrometric accuracy of the IRAM 
interferometer. Different values for T m from several transitions 
for the same molecule, or its isotopologues, may reveal vertical 
temperature gradients, or, in case of non-LTE excitation, den- 
sity gradients, since such transitions probe different regions of 
the disk. The values of T m from different molecules may also 
provide constraints on the density, because of the different crit- 
ical densities. And, more directly, the values of E,„, p m , R m and 
R out will reflect the chemical composition of the disk as func- 
tion of radius. 

4.2. Geometric and kinematic parameters 

Figure [7] shows the geometric and kinematic parameters PA, 
i and rotation velocity Vioo derived from the observed transi- 
tions in DM Tau, LkCa 15 and MWC 480. The grey area repre- 
sent the 1 <x range of the mean value computed from the CO 
lines. As expected, the derived value are in very good agree- 
ment altogether. In particular, the distribution of the derived 
values and the error bars are consistent with a statistical scat- 
tering, meaning that: i) the observed lines come from the same 
disk, ii) the error bars have the good magnitude, and iii) there 
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Fig. 7. From left to right: kinetic and geometric parameters derived for DM Tau, LkCa 15 and MWC 480. From top to bottom: 
rotation velocity at 100 AU, inclination and position angle. For DM Tau and MWC 480, the error bars represent (from left to 
right) the value of the parameter measured from the 12 CO J=2->1, 13 CO J=2->1, 13 CO J=1^0 and HCO + J=1^0 lines. For 
LkCa 15 (from left to right), the values are derived from the 12 CO J=2— > 1 line, from the simultaneous fit of the 13 CO J=2— > 1 and 
the 13 CO J=l— >0 lines, and from the HCO + J=l— >0 line. In all three sources, the grey area represent +lcr around the mean of 
the values of the CO isotopes. 



is no bias in the determination of the parameters. This analysis 
of the "basic" parameters is of prime interest, because it shows 
that the analysis method is robust, and that the error bars have a 
real physical meaning. Although noisier, the kinematic and ge- 
ometric parameters determined from HCO + are in agreement 
with those found from CO. We thus used the CO-derived val- 
ues to determine the surface density and excitation temperature 
ofHCO + . 



4.3. Results 

Tables [3] |4] and [5] present the best fit results for DM Tau, 
LkCa 15 a nd MWC 480, respectively. Compared to previous 
paper s (e.g. ISimon et alJ200ol: Duvertetal . 2000; Dutre v et al.l 
19981) . the new fits inclu de all the refinem ents of the method 
described in Sec|3]and in Pietu et al. ( 20061) . 

All quoted results where obtained using cylindrical co- 
ordinates and the continuum was subtracted to the spectro- 
scopic data before analysis. All the models were obtained in 
the same manner. In particular, we assumed a constant line 
width (e v - 0), Keplerian rotation (v = 0.5), and a flaring 
exponent e/, = - 1 .25 (the value for a Keplerian disk in hydro- 



static equilibrium with a temperature radi al dependence of 0.5). 
This v alue is very close to those found bv lChiang & Goldreich 
(119971) for the super heated layer model. Systematic tests with 
different sets of free parameters ensure that the geometrical and 
kinematical parameters are insensitive to these assumptions. 
The temperature T m (r) is also hardly affected. The surface den- 
sity S m (r) and outer radius R ollt depends to some extent on the 
assumed flaring parameter e/,. However, because of the cou- 
pling between scale height and line width, both Ho and dV 
should be taken with some caution. The dynamical mass de- 
rived from 12 CO also depends weakly on the assumed scale 
height, since we derive the inclinations from the apparent CO 
surface. 

The differences between the values repo r ted he re for 
DM Tau and those published by iDartois et al. ( 20031) com e 



from t he different assumptions in the analysis. IDartois et al 



(2003) assumed hydrostatic equilibrium, and model the line 
width following Eq|6] . They also assumed [ 12 CO]/[ 13 CO] = 
60. 

The 13 CO emission in MWC 480 is strong enough to al- 
low separate fitting of the J=2— >1 and J=l— »0 lines. We thus 
present separate results for each line. We also present a result in 
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Table 3. Best parameters for DM Tau. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Lines 


12 CO J=2^1 


13 CO J=2— >1 


13 CO J=l->0 


Mean 


HCO + I=l->0 


Systemic velocity, Vlsr (km.s~ ! 


) 6.026±0.002 


6.031 ± 0.003 


6.088 ± 0.004 


6.038±0.002 


6.01 ±0.01 


Orientation, PA (°) 


63.9 ±0.3 


65.7 ±0.4 


65.6 ±0.5 


64.7 ±0.2 


65 ± 1 


Inclination, i (°) 


-33.5 ±1.0 


-35.5 ±1.2 


-35.8 ±1.5 


-34.7 ±0.7 


-33 ±2 




Velocity law: V(r) 


" ^oo^iooAu) 








Velocity at 100 AU, V, 00 (km.s' 


- 1 ) 2.24 ±0.06 


2.07 ±0.06 


2.15 ±0.08 


2.16 ± 0.04 


2.10 ±0.15 


Stellar mass, M, ( M Q ) 


0.56 ±0.03 


0.48 ±0.03 


0.52 ± 0.04 


0.53±0.02 


0.5 ±0.1 


Surface Density at 300 AU, (cirT 2 ) 1 .4 ± 0.6 10 17 


5.3 ± 1.1 10 15 


5.9±0.610 15 




1.1 + 0.310" 


Exponent p 


3.8 ±0.3 


3.4 ±0.4 


2.8 ± 0.3 




2.3 ± 0.4 


Outer radius R ouU (AU) 


890 ±7 


720 ± 19 


760 ± 22 




800 ± 30 


Temperature at 100 AU, (K) 


26.0 ± 0.5 


14.5 ± 0.5 


8.0 ±0.3 




14 ±2 


Exponent q 


0.49 ± 0.01 


0.12 ±0.03 


-0.34 ± 0.07 




0.40 ±0.1 


dV (km.s- 1 ) 


0.12 ±0.01 


0.15 ±0.01 


0.20 ± 0.01 




0.13 ±0.01 


Scale Height at 100 AU, (AU) 


30+1.1 


28 ±4 


28 ±5 




(19 ± 8) 



Column (1) contains the parameter name. Columns (2,3,4) indicate the parameters derived from 12 CO J=2— >1, 13 CO 1=2— >1 and 13 CO J=l— >0 
respectively. Only 12 CO constrain the size of the inclination. Column (5) indicates the mean value of the kinematic and geometric parameters. 
Column (6) indicate the parameters derived from HCO + J=l— >0, obtained using a fixed inclination, orientation and stellar mass. For HCO + , 
because of the lower resolution, the error on the temperature and surface density include contributions from the uncertainty on the exponent. 



T m is better constrained at 300 AU (9.0 ± 0.5 K), and S m 
Table 4. Best parameters for LkCa 15 


at 500 AU (3.2 


+ 0.5 10 12 cnr 


- 2 ). 


(1) 


(2) 


(3) 


(4) 


(5) 


Lines 


12 CO L 


=2^1 


13 CO 


Mean 


HCO+ 1=1^0 


Systemic velocity, V LSR (km.s -1 ) 


6.29±0.01 


6.33 ± 


0.02 


6.30 ±0.01 


6.26 ± 0.04 


Orientation, PA (°) 


150.9 


±0.4 


150.4 


±1.0 


150.7 ±0.4 


150 ±2 


Inclination, i (°) 


51.8: 


t0.8 


49.6: 


tl.O 


51.5±0.7 


51.5 ± 1.3or46±3 


Velocity law: 


V(r): 


V 10f(l00Au) 








Velocity at 100 AU, V m (km.s- 1 ) 


3.07 4 


0.06 


2.97 + 


0.06 


3.00 ±0.05 


2.91 ±0.05 or 3.11 ±0.15 


Stellar mass, M, ( M Q ) 


1.06 ± 


0.04 


0.99 ± 


0.05 


1.01 ± 0.03 




Surface Density at 300 AU, (cirr 2 ) 


2.9 ± 0. 


1 10 16 


2.6 ± 0. 


210 15 




8.0 ± 0.5 10 12 


Exponent p 


4.4 ± 


0.3 


1.50 ± 


0.15 




2.5 ±0.3 


Outer radius R om , (AU) 


905 i 


z 40 


550 h 


= 20 




660 ± 60 


Temperature at 100 AU, (K) 


22.1 i 


z0.5 


21. 2d 


: 1.7 




[19] 


Exponent q 


0.37 ± 


0.02 


0.40 ± 


0.10 




[0.38] 


dV (km.s- 1 ) 


0.19 ± 


0.01 


0.29 ± 


0.03 




0.23 ± 0.04 


Scale Height at 100 AU, (AU) 


19 =1 


z 1 


17 d 


= 3 




(±) 



Column (1) contains the parameter name. Columns (2,3) indicate the parameters derived from 12 CO 1=2— >1, and a combined fit of the 
13 CO J=2— >1 and 13 CO J=l— >0 respectively. Column (4) indicates the mean value of the kinematic and geometric parameters. Column (5) 
indicates the HCO + results: we used a fixed temperature law compatible with the CO isotopologue results (see text). 



which both lines are fitted together: as expected in this case, the 
derived temperature is in between the results from the J=2— >1 
and J=l->0 transitions. However, for LkCa 15, the 13 CO emis- 
sion is too weak for a separate analysis: only the results from 
the simultaneous fitting of both 13 CO lines is presented. 

For HCO + in MWC 480, the signal to noise did not allow 
an independent derivation of the temperature T m : we used the 
law derived from 13 CO. A similar problem occurs in LkCa 15: 
for simplicity, we used T m = 19 K and q m = 0.38, val- 
ues compatible with the CO results, although there is a solu- 



tion marginally better (1.7<x) with slightly lower temperatures 
(T m = 13 K), with a somewhat steeper surface density gradient. 



5. Discussion 

We discuss here the physical implications of the results pre- 
sented in the preceding section. The following subsections fo- 
cus on the measurement of star masses, of the vertical tem- 
perature gradient, of the hydrostatic scale heights, of the CO 
abundances and CO outer radii, and of the HCO + distribution. 



Pietu et al.,: Vertical structure of protoplanetary disks surrounding HAe and TTauri stars 13 



Table 5. Best parameters for MWC 480. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Lines 


12 CO J 


=2^1 


13 CO J= 


=2^1 


13 CO J= 


=l->0 


13 CO 


HCO+ 


J=1^0 


Systemic velocity, Vlsr (km.s -1 ) 


5.076± 


0.003 


5.16 ± 


0.02 


5.17 ± 


0.02 


5.084 ± 


0.003 


5.11 


±0.04 


Orientation, PA (°) 


58.0 


±0.3 


55.8: 


t0.6 


57.1 : 


tO.9 


57.8- 


:0.2 


62 


±3 


Inclination, ( (°) 


37.5 


±0.7 


37.6: 


fcl.6 


33.2: 


t3.2 


37.0 


±0.6 


34.0 


± 1.5 




Velocity law 


V(r 


) - ^ioo( 100AU ) 












Velocity at 100 AU, V 100 (km.s- 1 ) 


3.95 i 


:0.06 


3.98 ± 


0.14 


4.45 + 


0.33 


4.03+ 


0.05 


3.88 


±0.13 


Stellar mass, M, ( M G ) 


1.76d 


:0.06 


1.65 ± 


0.15 


2.47±0.23 


1.83± 


0.05 


1.70 


±0.12 


Surface Density at 300 AU, (cm -2 ) 


3.6 ±0 


210 16 


3.3 ± 0. 


4 10 13 


5.3 + 


910 15 


4.5 ±0 


3 10' 5 


3.4±0.310 12 


Exponent p 


4.7 + 


0.3 


3.0 ± 


0.4 


4.5 ± 


0.9 


3.9 ± 


0.3 


1.5 


±0.2 


Outer radius R out , (AU) 


740: 


t 15 


450 + 


: 15 


520 i 


:70 


480- 


:20 


520 


±50 


Temperature at 100 AU, (K) 


48: 


t 1 


28 + 


:2 


21 i 


:4 


23 d 


: 1 


(15 


±6) 


Exponent q 


0.65 + 


0.02 


0.37 ± 


0.08 


0.28 + 


0.09 


0.37 ± 


0.04 


(0.6 


±0.4) 


dV(km.s-') 


0.25 + 


0.01 


0.25 + 


0.02 


0.15 + 


0.03 


0.21 ± 


0.02 


0.33 


±0.07 


Scale Height at 100 AU, (AU) 


10 ± 


1.1 


18 + 


-4 


20 i 


:2 


19 d 


:2 







Column (1) contains the parameter name. Columns (2,3,4) indicate the parameters derived from 12 CO J=2 — > 1 , 13 CO J=2 — > 1 and 13 CO J=l— >0 
respectively. Column (5) indicates the mean value of the kinematic and geometric parameters (Vlsr, PA, (', V,oo). For the other parameters, 
Column (5) indicate the results of a simultaneous fit of the two 13 CO transitions as for LkCa 15. Column (6) indicates the HCO + results: as the 
temperature is poorly constrained, the surface densities were derived assuming the temperature given by I3 CO J=l— »0. The surface density is 
best determined at 250 AU: 4.3 ± 0.3 10 12 cm 4 . 



5.1. CO dynamical masses 



Table 6. Stellar masses 



Table [6] shows a comparis on of the mass det ermination ob- 
tained in this work and in Simo n et al. (2000). The subtrac- 
tion of the continuum has slightly i mproved the mass deter- 
mination in the case of MWC 480. ISimon et aD d2000h used 
the old Solar metallicity (Z = 0.02) pre main-sequence (PMS) 
models. In FigJS] we present here also results for a metallicity 
Z = 0.01, much closer to the new Sola r metallicity (0.0126, 
see lAsplund et al.1 l2004t iGrevesse et al.l 120051 and references 



therein). LkCa 15 is in agreement with models of both metal- 
licities (M = 0.95 ±0.05 M G , t = 5-6 Myr and 1.1 ±0.1 M , t = 
4 - 6Myr, for Z = 0.01 and Z = 0.02 respectively). In fact, the 
measured dynamical mass does not provide any strong con- 
straint on the location of LkCa 15 in this diagram (since the 
PMS tracks are almost parallel to the Y axis for stars with 
0.5 < M»/M G < 1 in this distance independent HR diagram), 
but serves as a consistency check. However for MWC 480, ei- 
ther the metallicity is Z = 0.01 (i.e. nearly solar), in which case 
the de rived dynamical mass is in agreement with the lSiess et al. 
d200Ci tracks (M = 1.8M Q ,f = 8 Myr), or the metallicity is 
0.02 and MWC 48 would be lo cated at a larger distance as 
suggested by lSimon et al. (2000) (mass 1.9 - 2.0 M , age ~ 7 
Myr, and the distance would need to be increased by about 10 
%) 



5.2. Vertical Temperature Gradient 

The prim ary goal of these o bservations was to confirm the find- 
ings of Dartoi s et al.l d2003l) : the existence of a vertical temper- 
ature gradient. Fig|9]shows the temperatures deduced from the 
various CO lines versus the star effective temperatures in all 
sources of the sample. We confirm the DM Tau results with this 
slightly different analysis. MWC 480 also presents a significant 



Source 


Previous mass 


This work 




(Mb) 


(M Q ) 


DM Tau 


0.55 ± 0.03 


0.53 ± 0.02 


LkCa 15 


0.97 ± 0.03 


1.01 ±0.03 


MWC 480 


1.65 ±0.07 


1.83 ±0.05 



See Simon et al., 2000 for the previous measurements. The old masses 
are based only on 12 CO while the new measurements are the average 
of all available CO lines. 



temperature difference between 12 CO and 13 CO. However, in 
LkCa 15, there is no clear evidence of temperature gradient. 
As a general trend, the hotter stars of the sample, namely the 
Herbig Ae stars, exhibit larger kinetic temperature in the disks, 
not only at their surfaces but also for their interiors (see also 
Sect|5J]for AB Aur). For MWC 480, even if the surface ap- 
pears hotter, a significant fraction of the disk (for r > 200 AU) 
remains at temperature below the CO freeze-out temperature. 

The comparison between the three sources can be un- 
derstood by noting that the surface densities of 13 CO de- 
creases from DM Tau to MWC 480 and LkCa 15 (see Sect|5~8Tl. 
Moreover, the temperature is lower in DM Tau, resulting in 
higher opacities for the J=2— > 1 and J= 1 — >0 transitions than in 
the other sources. As a result, the location of the r = 1 sur- 
face in the 13 CO lines differ in the three sources. In DM Tau, 
this surface is above (for the J=2— >1 line) or at (for the J=l— >0 
line) the disk plane (as seen from the observer). A similar be- 
havior is observed in MWC 480. 

In LkCa 15, the 13 CO lines are nearly optically thin 
throughout the disk, and the temperature is accurately deter- 
mined by the J=2— »1/ J=l— >0 line ratio. The lack of differ- 
ence between the 12 CO and I3 CO results indicate both lines are 
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Fig. 8. Summary of the dynamical mass measurement for the sources of this sample. As in ISimon et al. (2000), the results are 
shown in a distance corrected HR diagram (L/M 2 versus T e ff)- New results have been obtained for MWC480 and LkCa 15. 
They are all compatible with the previous ones. In the case of MWC 480, the improvement in the mass determination is due to a 
better fit of the continuum, as explained in Sect l3.3.2l Le ft: Tracks for a ste llar metallicity Z = 0.01. Right: Tracks for a stellar 
metallicity Z = 0.02. The models have been computed bv lSiess et al. (2000). Dots starts at 1 Myr and are spaced by 1 Myr. The 
tracks range from 0.5 to 2.0 by 0. 1 M , with the 2.2 and 2.5 M Q tracks added. For T e ff, the error bars are assumed to be +/- one 
sub spectral type class. 



formed in a similar region: there must be little 13 CO at temper- 
atures much lower than indicated by 12 CO. 

5.3. CO Surface density and outer radius 

There are significant differences in surface densities between 
the three sources. At 300 AU, the measured surface densities 
of 13 CO range from ~ 2.5 10 15 ctrT 2 in LkCa 15 to ~ 3.2 10 15 
ctrr 2 in MWC 480 and ~ 7 10 15 ctrT 2 in DMTau. The varia- 
tions are even larger for 12 CO with values ~ 3 10 16 ctrT 2 in 
LkCa 15 and MWC 480, and 4 times larger, 14 10 16 ctrT 2 for 
DM Tau (although the latter is uncertain by a factor 2 because 
of the steep density gradient). 

The 12 CO/ 13 CO ratio is 20 ± 3 in DM Tau near 450 AU, 
the distance where this ratio is best determined. It is 1 1 + 2 at 
300 AU for LkCa 15, and 8 ± 2 for MWC 480. These values are 
much lower than the standard 12 C/ 13 C ratio in the solar neigh- 
bourhood. These low values are most likely the result of sig- 
nificant fractionation, as in the Taurus molecular clouds, since 
the typical temperatures near 400 AU are around 12 K (see also 
SectEl. 

The exponent of the surface density distribution varies sig- 
nificantly. For DMTau and MWC 480, an exponent of ~ 3.5 is 
appropriate for all isotopologues. For LkCa 15, there is a large 
difference between the exponent for 13 CO, 1.5, and that for 
12 CO, 4.5. Each line however samples a very different region 
of the disk. The 13 CO data are most sensitive to the 200 - 300 
AU region, while the 12 CO data is influenced by the outermost 
parts, 500 - 700 AU (because of the opacity of the 12 CO line) 
and up to the outer disk radius. The difference in exponent 
may thus indicate a steepening of the 12 CO surface density 
distribution at radii larger than > 500 AU. At such radii, the 
12 CO emission becomes optically thin, and better constrains 
the surface density than 13 CO. For DMTau and MWC 480, the 



surface density of CO falls faster with radius than usually as- 
sumed for molecular disks (p ~ 1 - 1.5, for H2). Note that 
the uncertainty on temperature law affects only weakly the sur- 
face density derivation, specially for the J=2— >1 transition for 
which the emission is essentially proportional to S m for the 
considered temperature range in the optically thin regime. This 
apparent steepening could reflect a change in the (H2) surface 
density distribution in the disk. However, it can also be a result 
of chemical effects. The isotropic ambient UV field will tend 
to photo-dissociate CO isotopologues near the disk outer edge. 
Another effect is that, because of the lower temperatures, de- 
pletion onto dust grains may be more effective at larger radii 
(although the lower densities result in a larger timescale for 
sticking onto grains). 

Table 7. HCO + and CO Outer Radius 



Source 


#„„,( 12 CO) 


fi ,„( 1J CO) 


HCO+ 




(AU) 


(AU) 


(AU) 


DMTau 


890 ±7 


740 ± 15 


800 ± 80 


LkCa 15 


905 ± 40 


550 ± 20 


660 ± 60 


MWC 480 


740 ± 15 


480 ± 20 


520 ± 50 



In all three sources, the outer radius in 13 CO is significantly 
smaller than in 12 CO. In LkCa 15, the large difference may be 
partly attributed to the change in the exponent of the density 
law. However, if we assume p = 1.5 for 12 CO, we get an outer 
radius of 710 + 10, still significantly larger than that for 13 CO. 
The effect is thus genuine. 



5.4. The LkCa 15 inner hole 

In the case of L kCa 15, it is also im portant to mention its pe- 
culiar geometry: iPietu et al. (2006) have discovered a 50 AU 
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Fig. 9. Temperatures derived from the CO isotopes versus ef- 
fective temperature of the central star. From left to right, the 
sources are DM Tau (filled triangles), GM Aur (empty stars), 
LkCal5 (filled squares), MWC480 (filled pentagons), HD 
34282 (empty stars) and AB Aur (filled hexagons). From top 
to bottom, 12 CO J=2->1, 13 CO J=2->1, 13 CO J=1^0 temper- 
atures. 



cavity in the continuum emission from LkCa 15. We checked 
whether there is any hint of this cavity in the CO data. Adding 
the inner radius as an additional free parameter yields Ri nt = 
13 + 5 AU for 12 CO and R int = 23 ± 8 AU for 13 CO. These 
values are consistent with a smaller hole, or no hole at all, in 
CO and I3 CO. The 50 AU inner radius determined from the 
continuum is excluded at the 7 cr level in 12 CO, and at the 3 <x 
level in I3 CO. Moreover the fit of a hole does not significantly 
change the value of p, for both 12 CO and I3 CO emissions. 

These small values indicate that the CO gas extends well 
into the continuum cavity. The cavity is thus not completely 
void, but still has a sufficient surface density for the I3 CO 
lines to be detected. The LkCa 15 situation is similar to that 
of AB Aur, for which the apparent inne r radius increases from 
12 CO to 13 CO and to the dust emission dPietu et al.ll2005l) . 



5.5. HCO + surface density 

The HCO + data analysis reveal 4 major facts 

- the outer radius in HCO + is similar to that in 12 CO, and in 
particular it is larger than that in 13 CO. 

- in the two cases where it could be determined, the excita- 
tion temperature of the HCO + J=l— >0 is lower than that of 
the I3 CO J=2^1 line. 

- at 300 AU, the [ 13 CO]/[HCO + ] ratio is 500±150 (DM Tau), 
320 ± 50 (LkCa 15), and 1300 ± 200 for MWC480, with 
an additional uncertainty around 50 % for these last two 
sources due to the adopted temperature distribution. These 
values ar e similar to, although sli ghtly smaller than the one 
found bv lGuilloteau et al" ( 19991) for the circumbinary disk 



- This ratio varies with radius: there is a tendency to have a 
somewhat flatter distribution of HCO + than of CO (with 
the exception of 13 CO in LkCa 15). 

In DM Tau, the excitation temperature falls down to about 6 K 
beyond 600 AU, while the CO data indicates kinetic temper- 
ature of about 10 K in this region. Our calibration technique 
ensures that this is not an artifact. It may be a hint of very low 
temperatures in the disk mid-plane. On the other end, it may be 
an indication of sub-thermal excitation. However, as the criti- 
cal density for the HCO + J=l— >0 is only a few 10 4 cm -3 , such 
a sub-thermal excitation would require that the HCO + gas is 
located well above the disk plane. Spatially resolved observa- 
tions of the HCO + J=3— »2 transition are required to check this 
possibility. 

5.6. Turbulence in Outer Disks 

We derive intrinsic (local) line widths ranging between 0.12 
and 0.29 km.s . When taking into account the thermal compo- 
nent (0.08 to 0.15 km.s" 1 ), from Eqs|6]and|7] we derive turbu- 
lent widths below 0.15 km.s 1 . These values should be used as 
upper limits, since the spectral resolution used for the analysis 
(0.2 km.s -1 ) is comparable to the derived line widths. They 
are nevertheless significantly smaller than the sound speed, 
C s = 0.3 to 0.5km.s _1 in the relevant temperature and radius 
range. The turbulence is thus largely subsonic. A more pre- 
cise analysis, using the full spectral resolution and an accurate 
knowledge of the kinetic temperature distribution, is required 
for a better determination. 

5.7. Hydrostatic Scale Height 



(2006) 



of GG Tau (~ 1500). On the other end, iPetv et al. 
find a ratio of about 10 4 for the low mass disk around 
HH 30, but at a smaller radius (100 AU). 



We follow the approach which is described in Sect l3.2.5l where 
the scale-height is not formally derived from the hydrostatic 
equilibrium. In a first step, the flaring exponent <?/, is taken 
equal to e/, = -1.25. Such a value is in agreement with pre- 
dictions from models s uch as lChiang & Goldreichl dl997) and 
D' Alessio et al. d 1999b . The observed transitions indicate ap- 
parent scale heights at 100 AU of 30 AU in DM Tau, and 19 AU 
in LkCa 15 and MWC480 (see Tab EE). The only discrepant 
result is the value derived from 12 CO in MWC480, 10 AU. In 
a second step, treating <?/, as a free parameter provides similar 
results, except for LkCa 15 where the I3 CO data is best fitted 
by a flat (e/, = 0), but thick disk (see Fig|5]l. 

On the other end, using the temperature derived from the 
13 CO J=2— >1 transition as appropriate, we can derive the hy- 
drostatic scales heights. These scale heights are 15, 13 and 9.5 
AU, respectively for DM Tau, LkCa 15 and MWC480, with 
a typical exponent e/, = -1.35, approximately 1.5-2 times 
smaller than the apparent values. Although, as mentioned in 
Sect l3.2.5l the scale height has to be interpreted with some 
caution (especially when determined from I2 CO), these results 
suggest that CO has a broader vertical distribution than what is 
expected from an hydrostatic distribution of the gas. The cur- 
rent data do not allow us to identify the causes of this effect. 
Vertical spreading due to turbulence can be rejected as the tur- 
bulence is largely subsonic (see Sect l5.61 >. The larger apparent 
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height may be due to CO being more abundant in the photo- 
dissociation layer above the disk plane. 

5.8. Comparison with other Disks and Models 

Only one other source has been stu died at the same level 
in the CO isotopologues: AB Aur by [Pietu et all d2005l) . The 
temperature behavior of the ABAur disk is similar to that 
of the MWC 480 disk, although the AB Aur disk is warmer. 
Temperature derived from CO for all sources are summarized 
in Figj9] As expected, the disk temperatures (midplane and sur- 
face) clearly increase with the stellar l u minos ity. The A0/A1 
star HD34282, mapped by IPietu et al.l d2003l) in 12 CO also 
exhibits a hotter disk surface while the TTauri star GM 



Aur (Dutrey et al. 1998) has a colder disk. From Table 1 of 
Piet u et al.1 d2005l) . the 13 CO surface density in AB Aur at 300 



AU is about 210 cm , 4 to 8 times larger than in the DM Tau, 
MWC 480 or LkCa 15 disks. This larger surface density prob- 
ably results from a lack of depletion onto dust grains, since 
the temperature in the AB Aur disk remains well above the CO 
condensation temperature at least up to 700 AU in radius, con- 
trary to those of DM Tau, LkCa 15 and MWC 480 which are 
much colder beyond 200 - 300 AU. 

LkCa 15 differs slightly from the general picture: no clear 
vertical temperature gradient is observed, and the 13 CO surface 
density slope is much less steep (exponent p = 1 .5) and the 
CO isotopologues indicate a discrepant scale height behaviour. 



This m ay be linked to a peculiar disk geometry: IPietu et al 



(2006) tave reported a 50 AU cavity in the continuum emis 
sion from LkCa 15, perhaps due to planetary formation or a 
low mass companion. This peculiar geometry may affect the 
temperature structure of the surrounding disk, by for example 
shielding it from stellar radiation behind a warm, thick, inner 
rim. Any definite conclusion about the peculiarities of LkCa 15 
seems premature at this point. 

Molecular abundances in proto-planetary disks have been 
studied by several authors. Disks are believed to be dominated 
by a warm molecular layer, exposed to the UV radiat ion ema- 
nating from the central star. Ivan Zadelhoff et al.l (12003b showed 
2-D radiative transfer effects result in a stronger penetration of 
the UV field into the disk than simple 1+1 D models assumed. 
A mor e refined model was published by lAikawa & Nomura 
(2006), which use a time dependent chemistry taking into ac- 
count the thermal structure of the disk, 2-D radiative trans- 
fer and sticking onto grains. Their model parameters are more 
appropriate for DM Tau. Despite its comprehensiveness, this 
model fails to represent the observed surface densities: it pre- 
dicts relatively shallow (p ^ 1) radial distribution of the sur- 
face density, in sharp contrast with the steep slope observed. 
However, their model only extends out to 300 AU, while the 
constraint we obtain on the surface densities in DM Tau refer 
mostly to larger radii. 



Semenov et al .1 (120061) point out the importance of vertical 



mixing in explaining the mo lecular surface den sities. A simi- 
lar study was conducted by lWillacy et al. (2006) who conclude 
that a diffusion coefficient of order 10 18 cm 2 s _1 is required 
to bring the CO surface density in reasonable agreement with 



observations. Vertical mixing may be essential in explaining 
why CO is detected at temperature well below the freeze out 
value, and also why, despite very different temperature profiles, 
DM Tau and MWC 480 display similar CO surface densities. 
None of these model studied the 13 CO isotopologue. 



Hily-Blant et al. (2007) used a simpler PDR model to study the 
effects of grain size and stellar UV flux on the CO isotopologue 
abundances. In this static model, self-shielding of H2, CO and 
13 CO is treated explicitly, but CO sticking onto grains is not in- 
cluded. This model shows the importance of fractionation and 
reproduces better the observed surface densities. In particular, 
from their Fig.6, [ 12 CO]/[ 13 CO] ratio of order 10 - 20 is ob- 
tained for large grains (a+ = 3 10~ 4 m) and^- = 1000 (x being 
the enhancement over the interstellar UV field at 100 AU, with 
a 1/r 2 dependence), i.e. when the UV field penetration is im- 
portant. 

So far, no published model reproduces the smaller outer 
radii for 13 CO than for 12 CO. To do so, it seems essential to ac- 
count for selective photodissociation of CO isotopologues by 
the interstellar UV radiation field, and in particular by the UV 
radiation impinging on the disk outer edge. The disappearance 
of 13 CO at distances ~ 500 AU from MWC 480 and LkCa 15 
indicate that the typical extinction through the disk at this ra- 
dius corresponds roughly to Ay - 1, although the disks extend 
at least to 700 - 900 AU. A consistency check can be done 
on this hypothesis: assuming a standard Ay/Hi ratio, the 13 CO 
outer radius can be used to constrain the CO abundance at this 
distance. From Fig|U|6] the typical 12 CO surface density at the 
13 CO outer radius is 2-4 10 15 crrT 2 , which converts to X[ 12 CO] 
^3-6 10~ 5 , in good agreement with the expectation in this 
transition region. In DM Tau, Ay — 1 would occur near 800 
AU, i.e. the (outer) DM Tau disk is more opaque than those of 
MWC 480 and LkCa 15. This larger opacity can result either 
from a larger H2 density, or from a larger Ay/Hi, i.e. more 
small grains in DM Tau. 

All the chemical models are hampered by the need to spec- 
ify the appropriate H2 density distribution and the dust proper- 
ties. These parameters remain largely unknown. The values of 
p and R nnt for the dust in LkCa 15 and MWC 480 derived by 
Pietu et all d2006l) differ systematically from those found here 
for CO. In addition, the absolute value of the absorption co- 
efficient of dust is also uncertain, which makes the derivation 
of an H2 surface density through contin uum observations q uite 
difficult. Using the parameters given by lPietu et al.l (|2006), we 
find at the dust core outer radius, ^ 180 AU, that the CO abun- 
dance is ~ 10~ 6 for LkCa 15, while it is ~ 3 10~ 7 -3 10~ 6 for 
MWC 480, depending on which temperature profile is adopted 
for the dust. These very low values, as well as the very differ- 
ent values of the exponent p for dust and CO, may thus indicate 
that the H2 surface density cannot be extrapolated by a simple 
power law, but displays instead a very significant change near 
~ 150 AU. This can have significant impact on the chemical 
model predictions for any specific disk. 

6. Summary 

We report new 13 CO observations in LkCa 15 and MWC 480 
obtained with the IRAM array, and HCO + observations of 
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LkCa 15 and MWC 480 and DM Tau. These data were analyzed 
together with existing 12 CO data. As we improved our mod- 
elling method, we also carried again the analysis on DM Tau 
data as a consistency check. By comparing these results with 
previous one obtained on AB Aur (and to the temperatures de- 
rived from 12 CO observations in HD34282 and GM Aur), this 
allowed us to derive the following conclusions: 

- We systematically detect radial gas temperature gradi- 
ents. Vertical temperature gradients are also ubiquitous in 
ABAur, DM Tau, and MWC 480, but LkCa 15 is an ap- 
parent exception (perhaps due to its peculiar geometry). 
Disks surrounding hotter stars, such as Herbig Ae objects 
(MWC 480 and AB Aur), present a hotter CO surface. This 
behavior agrees with the prediction from models of passive 
irradiated disks. 

- CO disks around hotter stars are also hotter around the 
mid-plane (ABAur and to a lesser extent MWC 480). 
Although [CO/H2] ratios are not directly measurable, the 
much higher surface density of 13 CO in the warmest disk, 
AB Aur, indicates that depletion onto grains is an important 
mechanism for the colder disks. However, the CO surface 
density does not appear to be a monotonous function of the 
stellar type. 

- For the more sensitive transitions, we are able to measure 
an apparent scale height which appears to be larger than 
the hydrostatic scale height. Although marginal, this result 
is in agreement with the expectation that molecules prefer- 
entially form in a warm photo-dissociation layer above the 
disk plane. 

- However, in the outer part of TTauri disks (R > 300 AU) 
and even around MWC 480, the averaged kinetic temper- 
ature measured along the line of sight is well below the 
freeze-out temperature of CO, but a significant amount of 
CO is still observed. This is in contradiction with most pre- 
dictions of chemical models but would be in favor of mod- 
els with vertical mixing. 

- The CO surface density falls off much more rapidly with 
radius than predicted by most chemical models. In the cases 
of DM Tau and LkCa 15, we also observe that the 12 CO has 
a steeper surface density exponent than 13 CO. This suggests 
a steepening of the surface density slope beyond 300 - 500 
AU, but photo-dissociation, and also perhaps depletion on 
grains are likely to play a role. 

- A change in surface density behaviour is also suggested by 
the very different parameters (exponent p and outer radius) 
derived from the continuum emission. 

- For all disks where 12 CO and 13 CO data are available, 
photo-dissociation by the ambient UV field impinging on 
the disk edge is likely the appropriate mechanism to ex- 
plain the smaller outer radii determined in 13 CO compared 
to that of 12 CO. 

- Except for the warmest disks, CO fractionation is an im- 
portant mechanism. For DM Tau and LkCa 15, we estimate 
an isotopologue ratio, of order or less than 10-20 between 
300 AU and the 13 CO outer radius. 

- At 300 AU, the [ 13 CO]/[HCO + ] ratio is around 600, within 
a factor 2. It apparently decreases with radius. 
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